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We investigate the mechanism of the resonant X-ray scattering on the Ti K edge region of YTiOa using the band structure 
calculation combined with the local density approximation. A large intensity is obtained for the orbital superlattice spots. The 
calculated spectra consist of several peaks as a function of photon energy in agreement with the recent experiments. Against 
a naive interpretation that directly relates the intensity to the orbitally polarized 3d states, the obtained large intensity arises 
from the distorted crystal structure, i.e., the tilt of the TiOe octahedra and the Jahn- Teller distortion, which considerably 
modifies the 4p states in the intermediate states of the dipolar process. This casts doubt on a prevailing assertion that the 
resonant x-ray scattering is a direct observation of the orbital order. 
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I. INTRODUCTION 

Recently the resonant x-ray scattering spectroscopy 
has attracted much attention because this experimen- 
tal method is considered as one of the most powerful 
tools which can directly observe the orbital order in the 
strongly correlated systems. A lot of experiments have al- 
ready been <[^ripf^9trt the perovskite transition-metal 
compounds,y'&m'§Ll V2O3EI and DyBaCa,! in which the 
scattering intensities were observed on the orbital super- 
lattice spots, which are forbidden for usual x-ray scat- 
tering. Such anomalous x-ray scattering will simply be 
referred to as RXS. Despite much effort so far, how and 
to what extent the RXS spectra are reflecting the orbital 
order are still controversial, especially on the perovskite 
compounds, since the orbital ordei- usually accompanies 
the Jahn-Teller distortion (JTD)E3 and the GdFeOa type 
distortion, which can also gives rise to the RXS intensity. 

The resonant scattering may be described by a sec- 
ond order process that a photon is virtually absorbed by 
exciting a core electron to unoccupied states and then 
emitted by recombining the excited electron with the 
core hole. Therefore, for the K edge of the transition- 
metal compounds, the 4p states are involved in the in- 
termediate states of the dipolar process. For the RXS in 
the K edge of LaMnOa, the so called Coulomb mecha- 
nism was proposed that the anisotropic part of the intra- 
atomic Coulomb interaction between the orbitally po- 
larized 3d states and the 4p states is responsible to tha 
RXS intensity by making the degenerate 4p levels split .til 
However, subsequent studies based on the band struc- 
ture calculations have revealed that the Coulomb effect 
is much smallef, [(4^/100) than the JTD effect on the 
RXS intensity.EJOllJ This is because the ip states are so 
extended in space that they are sensitive to neighboring 
O potentials. 



The RXS spectra on YVO3 have been observed for the 
pre K edge and the main K edge regions of V.H The 
spectra consist of several peaks as a function of photon 
energy, and the shapes are different between the (100) 
and (Oil) reflections. Assuming that the scattering ten- 
sor is directly related to the orbital polarization of the 3d 
states, the authors argued that the azimuthal angle de- 
pendence of ±he spectra was consistent with the C-type 
orbital ordeulj in the low temperature phase. However, 
this assumption made there implies the dominance of the 
Coulomb effect, which is questionable. Recently the res- 
onant x-ray scattering experiment has also been carried 
out on YTiOa, and the RXS spectra have been observed 
with their shapes similar to those of YVO3.I3 In this pa- 
per, concentrating our attention on YTiOs, we study the 
mechanism of the RXS spectra by carrying out a band 
structure calculation on the experimentally determined 
crystal structure.E£l 
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FIG. 1. Sketch ofthe crystal structure of YTiOs. The large 
spheres represent the oxygen ions, the small solid spheres the 
yttrium ions, and the octahedra the titanium ions. The ar- 
rows assigned to each Ti ion define the local ^ryc^ axes pointing 
to nearest neighbor oxygen sites. 

The crystal structure of YTiOs is known to belong 
to the space group Pbnm, as shown in Fig. 0; four Ti 
ions are inequivalent in the unit cell. The Ti-O-Ti bond 
angles strongly deviate from 180° due to the tilt of the 
TiOg octahedra,0 and Ti-0 bond lengths are different in 
TiOg octahedra due to the JTD. One of the t2g levels is 
occupied under the crystal field in the d} configuration of 
Ti3+ 

ions, and since the t2g states are weakly hybridizing 
with O 2p states, the JTD is smaller (~ 20%) than that 
on LaMnOall^. Note that the crystal structure of YVO3 
belongs to the same space group of Pbnm wit h || Si piilar 
sizes of the tilt of VOe octahedra and the JTDBI13. 

We calculate the RXS spectra using the band calcu- 
lation based on the local density approximation (LDA) 
together with the muffin-tin approximation. This ap- 
proximate scheme leads to a metallic ground state for 
YTiOa with a very small orbital polarization. In spite of 
this shortcoming for the 3d bands, we expect that the Ap 
bands are well described in the present calculation, since 
they have energies ~ 15 eV higher than the 3d bands 
and thereby the details of the 3d bands are irrelevant. 
What is important, which will become clear, to make the 
4p states different at different Ti sites is the O potential 
of neighboring sites. The Ap states are so extended in 
space that they are sensitively modified by hybridizing 
and interacting with the electronic states of neighbor- 
ing ions. Actually, for LaMnOa, it is known that the 
LDA+C/ method, which predicts an insulating ground 
state with large orbital polarization, gives the 4p bands 
and the RXS spectra nearly the same as those given by 
the LDA method, although the latter predicts a small 



band gap and a small orbital polarizationB 

We obtain the RXS spectra with sufficient intensities 
for YTiOs. The obtained spectra as a function of photon 
energy resemble closely to the results of the experiment 
for the (100) and (Oil) reflections of YVO3J and also 
to those of the recent experiment of YTiOsfl Since the 
muffin-tin approximation averages the Coulomb interac- 
tion between the Ap and the 3d states, the anisotropy of 
the Coulomb interaction is eliminated in the calculation 
of the RXS spectra. Therefore the obtained RXS inten- 
sity arises from the distorted crystal structure. We may 
roughly estimate that the Coulomb effect is less than 1 /4 
of the JTD effect, since the Coulomb effect is evaluated as 
1/100 smaller than the JTD effect in LaMnOa.BEa Note 
that the JTD effect cannot be distinguished from that 
of the tilt of the TiOg octahedra in the present calcula- 
tion. Such sensitivity of the 4p states to the electronic 
states at neighboring sites has already been recognized 
in the analysis of the magnetic circular dichroism of the 
K edge absorpticp. spectra on the ferromagnetic states of 
transition metalsEj and the analysis of. the resonant x-ray 
magnetic scattering spectra on CoOE3 and NiO.a 

This paper is organized as follows. In the next section, 
we briefly discuss the band structure calculation on the 
ground state. In section 3 we outline the calculation pro- 
cedure of RXS spectra and discuss the calculated results. 
The final section is devoted to concluding remarks. 
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CALCULATION OF THE GROUND STATE 
AND DENSITY OF STATES 
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FIG. 2. The 4p wave functions of an atomic Ti'^"'' ion, calcu- 
lated within the Hartree-Fock approximation. Solid line and 
dashed line represent the radial wave functions rR{r) with- 
out and with the core-hole potential. The arrows indicate 
the mufhn-tin radius tmt of Ti site and the typical distance 
dTi-o between Ti nucleus and O nucleus in YTiOa. 



First we examine the meaning of the atomic 4p states. 
Figure |^ shows the atomic 4p wave function for a Ti'ri| 
ion, calculated within the Hartree-Fock approximation.c3 
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The 4p wave function extends considerably in space so 
that a large part of the amplitude would lie on the neigh- 
boring oxygen sites or the interstitial region in solids. 
This situation is not changed in the presence of the core- 
hole potential, as shown by the broken line in the figure. 
Thus one may easily understand the sensitiveness of the 
Ap states to electronic states at neighboring sites. The 
concept of the atomic 4p levels loses its clear meaning 
in solids, and the ip states are to be treated as energy 
bands. Nevertheless, the p symmetric states (centered at 
a Ti site) above the Fermi energy are referred to the Ap 
states, since the atomic character is almost preserved in- 
side the atomic sphere in the band structure calculation. 

In the band structure calculation for YTiOs, we use 
the Korringa-Kohn-Rostoker (KKR) method within the 
LDA scheme. We assume the crystal structure detet- 
mined experimentally at room temperature in Ref.EEl, 
which is shown in Fig. |l|. As mentioned before^-Jthe crys- 
tal structure belongs to the space group Pbnm-jt3 the unit 
cell involves four Y ions, four Ti ions and 12 O ions. The 
four Ti sites are distinctive in terms of the orientation 
even though the orbital ordering does not occur. The 
large tilt of TiOg octahedra makes the Ti-O-Ti bond an- 
gles strongly deviates from 180° by as much as 35°. Each 
TiOe octahedron also undergoes the JTD such that the 
direction ^ corresponds to the short bond, the rj to the 
long bond, and the C to the middle bond, where three 
arrows ^77^ attached to each Ti ion are directed to neigh- 
boring oxygen sites and are almost orthogonal to each 
other. There exists an important symmetry that the con- 
figuration surrounding the Ti 1 site becomes coincident 
to that surrounding the Ti 2 site under the rotation of 
all positions of atoms by tt around the line parallel to the 
X axis passing through the Ti 1 site. In the same way, 
the configuration surrounding the Ti 1 site becomes co- 
incident to those surrounding the Ti 3 and the Ti 4 sites 
under the rotation of the atomic position by tt around the 
line parallel to the y and the z axes, respectively, passing 
through the Ti 1 site. 

As regards the magnetism, the system undergoes a 
ferromagnetic transition below the Curie temperature 
Tq — 29K.Ej Neglecting the spin-orbit couphng, we ob- 
tain a ferromagnetic ground state with the magnetic mo- 
ment 0.95/iB per Ti'^+ site (0.68/xb inside the muffin- 
tin sphere n£.jlj'^+ site), consistent with the previous 
calculations .t§t3 This value|4s a little larger than the 
experimental value 0.84yUB The orbital polarization 
is underestimated to be very small. The ground state 
is obtained as a metal, although it is experimentally a 
semiconductor. The LDA is known to have a tendency 
of wrongly predicting a metallic ground state. We think 
these shortcomings for the 3d states have only minor ef- 
fects on the ip states, since they have high energy (~ 15 
eV above the Fermi level). 
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FIG. 3. The DOS projected onto the s, p and d symmetric 
states inside the mufhn-tin sphere of a Ti site. The origin of 
energy is at the Fermi level. The d DOS is overwhelming the 
s and p DOS's. The inset shows the magnified DOS projected 
onto the p symmetric states. 

Figure ^ shows the density of states projected onto the 
s, p, and d symmetric states inside the muffin-tin sphere 
of a Ti site. Around 7 eV below the Fermi level, the 
3d states form bonding states with the 2p states of the 
neighboring O ions. Around the Fermi level, the 3d states 
are highly concentrating in a form of antibonding states 
with the O 2p states. In high energy region, a broad 
band appears, to which the s, p and d symmetric states 
equally contribute. The p DOS in this energy region is 
much smaller than the d DOS near the Fermi level, as ex- 
pected from the character of the atomic 4p wave function 
that extends considerably in the interstitial region. The 
p DOS is enlarged in the inset. It has two peaks. Note 
that the p DOS is roughly proportional to the K edge ab- 
sorption spectrum, if the effect of the core hole potential 
is neglected in the final states. This final state interac- 
tion may enhance the peak of the low energy side in the 
p DOS. In any event, the two-peak structufe corresponds 
well to the absorption spectrum of YVO3D 



III. CALCULATION OF RXS SPECTRA 
A. Formulation 

In the resonant process around the Ti K edge region, 
the Is core electron is virtually excited to the Ap states. 
Therefore the RXS spectra directly reflect the Ap states at 
the scattering center. We specifically calculate the RXS 
intensities for the (100), (001) and (Oil) reflections. 

We choose the x, y and z axes such that they are par- 
allel to the crystal a, b, and c axes, respectively. With 
respect to these axes, we define the Is-Ap dipole transi- 
tion density matrix rmn(£) at the Ti k site in the unit 
cell as 
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[i?L(r'KK0,.fe(r')]'^(£-£,fc), (3.1) 



b,k 



where 0^ represents the wave function with the band 
index b, wave- vector k and energy ^ . The Vm is a pro- 
jection operator which projects the wave function (j)^ 
onto the m (m G x,y,z) component of the p symmetric 
part around the Ti k site, and Ris represents the Ti Is 
wave function. This matrix js obtained from the KKR 
Green's function gmn{r, r'; £)E3 by the equation 
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(3.5) 



Using these representations, the scattering amphtudes 
are expressed in simple forms. For the (100) reflection, 
we have 



a{e) 
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(3.6) 



(3.2) 

Numericafly this ls-4p transition density matrix is 
roughly proportional to the 4p density matrix pmn(e) — 
— ilm J^'^'^ r'^dr gmnir,r; e) with vmt being the mufhn- 
tin radius. Assuming that photons are linearly polarized, 
we have the scattering amplitude at the Ti k site as 
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where A is a numerical constant, E^^ and -B™* are the 
771 and n components of polarization of the incident and 
the emitted photons, respectively. In the denominator, 
Lu is the photon energy, and Sis is the energy of the Is 
state. The F describes the broadening due to the core- 
hole lifetime, which are assumed 7^ = 1 eV. The effect 
of the core-hole potential is neglected in the intermedi- 
ate states. We think this effect on the RXS intensity is 
small, since the core hole potential is spherically sym- 
metric, and it affects merely the diagonal elements of the 
density matrix which are irrelevant to the RXS intensity. 
The integration in terms of e is performed from the Fermi 
energy to 60 eV above it. 

It is obvious that the transition density matrix 

(k) 

Tmn becomes rSymmetric with neglecting the spin-orbit 
interaction.c3'c3 Using a, b, c, a, f3 and 7, functions of e, 
the matrix at the Ti 1 site can be represented as 



.(1) 




(3.4) 



Due to the crystal symmetry mentioned before, the tran- 
sition density matrices r'^', t^^^ and t''^^ at the Ti 2, Ti 
3, and Ti 4 sites are represented as 



.(2) _ 




In the similar way, we have the scattering amplitude for 
the (001) and (Oil) reflections as 
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(3.7) 



(3.8) 



Only one of the off-diagonal elements remains in each ex- 
pression of the scattering amplitudes, with the diagonal 
elements a, &, and c completely vanishing. Therefore, we 
can select out one of them by choosing the reflection in- 
dex. The photon energy dependenc e is factored out from 
the polarization of photon in Eqs. ( ^.6D , (3.7), and (^). 
The RXS intensities /, 



(100) J 



1/(100)1 I 1/(001)1 7 and |/(oi 



^(001), and /(on) are given by 
) P, respectively. 



B. Calculated results 

We calculate the RXS spectra on the crystal struc- 
ture determined at room temperature in Ref.lZl, using the 
above formulas. We show in the following the results 
calculated in the ferromagnetic ground state of YTiOs, 
although the RXS experimental data are given above the 
Curie temperature Tc.Q We checked that the RXS spectra 
are insensitive to the magnetic states by calculating the 
spectra also in a nonmagnetic state. This is consistent 
with the 4p DOS above 15 eV in Fig. ^ which shows no 
exchange splitting. 
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FIG. 4. RXS spectra as a function of photon energy. The 
left panels are the calculated spectra for the (100), (001), and 
(Oil) reflections on YTiOa. The origin of energy corresponds 
to the photon energy of exciting an electron from the Is state 
to the Fermi level. Right panels are the experimental spectra 
in the low temperature phase of YVO3 (ref. 6) . 

Figure ^ shows the calculated spectra as a function 
of photon energy for the (100), (001) and (Oil) reflec- 
tions. Note that the anisotropy of the Coulomb interac- 
tion between the 4p and the 3d states is eliminated in 
the calculation by using the muffin-tin approximation. 
The experimental spectra in the low temperaturej-phase 
of YVO3 are also shown for the sake of comparison Jj since 
their shapes are quite similar to the unpublished experi- 
mental data for the RXS spectra of YTiOs.El As already 
mentioned, the resemblance of the RXS spectra is quite 
understandable in the light of the mechanism of the dis- 
torted crystal structure, since both the JTD and the tilt 
of the VOg octahedra in the low temperature phase of 
YVO3 are quite similar to those in YTiOa. In the (100) 
and (001) reflections, the intensities vanish for the <t <t' 
channel. For the cr — > tt' channel, we have three peaks 
denoted by A, B, and C, whose intensities are comparable 
to each other in the (100) reflection, while we have two 
peaks denoted as D and E in the (001) reflection. In the 
experimental data of YTiOa, three peaks appear at 4.974 
keV, 4.986 keV, and 5.0 keV, corresponding to peaks A, 
B, C, for the (100) reflection, while two peaks appear at 
4.974 keV and 4.984 keV with the former intensity about 
twice as the latter, corresponding to peaks D and E.l1 In 
the (Oil) reflection, both the a a' and the a —* tt' 
channels contribute to the intensities; we have a domi- 
nant peak denoted as F and a small peak denoted as G. 
The energy of peak F is found slightly higher (^ 2 eV) 
than that of the peaks A and D. In the experimental data 
of YTiOa, two peaks appear at 4.976 keV (2 eV higher 
than the positions of peaks A and D) and 4.986 keV with 
the former intensity dominating the latter. 



FIG. 5. Azimuthal angle dependence of the RXS intensities 
for the (100), (001) and (Oil) reflections. 

Figure]^ shows the azimuthal an gle de pen dence o f th e 
RXS spectra. As clear from Eqs. ( |3.6| ), (3.7), and (^^), 
the peaks belonging to the same scattering vectors show 
the same dependence. In the (100) and (001) reflections, 
the intensities are proportional to cos^ for the a — > tt' 
channel, with the condition E™ || & at '0 = 0. In the the 
(Oil) reflection, the intensities are proportional to sin^ 2'0 
for the cr — > cr' channel, with the condition E™ || a at 
"0 = 0, while the dependence is rather complicated for the 
CT — > tt' channel. The calculated results are-jn agreement 
with the those of YVOsB (and also YTiOaQ). 
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FIG. 6. The 4p density matrix pmn(e) at the Ti 1 site, as 
a function of energy e, (a) in the xyz coordinate system, and 
(b) in the ^t^C^ coordinate system. The thick lines represent 
the off-diagonal elements, and the thin lines represent the 
diagonal elements. The origin of energy is at the Fermi level. 

To see more closely how the Ap states are modified, we 
examine the 4p density matrix pmn(£) at the Ti 1 site, 
which is shown in Fig. |6[ The left panel shows the rep- 
resentation in the xyz system. Although the diagonal 
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elements are larger than the off-diagonal elements, the 
former quantity does not contribute to the RXS inten- 
sity. The right panel shows the representation in the local 
^Tyf^ coordinate system. The off-diagonal elements remain 
comparable to those in the xyz system, indicating that 
the ^?7C axes are still deviating from the principal axes of 
the density matrix. Since the off-diagonal elements are 
fluctuating with changing energy, the principal axes may 
also be fluctuating. One may think that the principal 
axes of the density matrix are related to the direction of 
the 2>d orbital polarization. However, this is not the case 
in the present calculation, since the Zd orbital polariza- 
tion has little effect on the 4p states within the muffin-tin 
approximation. We interpret this as the sensitiveness of 
the Ap states to the tilt of the TiOe octahedra which 
makes the local S^rjQ coordinate axes no longer principal 
axes. 



IV. CONCLUDING REMARKS 

We have studied the effect of the crystal structure to 
the RXS spectra of YTiOa, using the band structure 
calculation within the LDA. The calculated spectra are 
mainly compared with the experimental data of YVO3, 
since the unpublished data of YTiOa are similar to those 
of YVO3. This similarity is consistent with the mech- 
anism of the distorted crystal structure. The obtained 
spectra as a function of photon energy are in agreement 
with the experiments of YVO3 (and also YTiOa). Since 
the muffin-tin approximation averages the Coulomb in- 
teraction between the Ap states and the 3d states, the 
anisotropy of the Coulomb interaction is eliminated in 
the calculation of the RXS spectra. We have clearly 
shown that along with the JTD the tilts of the transition- 
metal oxide octahedra need to be taken into account. 
The agreement between the calculation and the experi- 
ment suggests that the RXS intensity originates mainly 
from the distorted crystal structure, casting doubt on a 
prevailing assertion that the RXS is a direct observation 
of the orbital order. We hope that the present work could 
accelerate further study on the relation between the RXS 
and the orbital order. 

In spite of the above findings, there are several points 
to be made clear. First, since the crystal structure in- 
cludes both the JTD and the tilt of the TiOe octahedra, 
it is not clear which effect is important. Closely related 
to this question, the crystal structure of YVO3 changes 
into a state of a very small JTD for T > 100 K (the 
tilt is not changed so much). According to this change, 
the intensities of peaks a and c were found to become 
very small, while those of other peaks were found to re- 
main similar. This suggests that the effect of the JTD 
on the RXS spectra is different from that of the tilt of 
the octahedra. For clarifying this point and for a direct 
comparison with the experiment, the calculation of the 
RXS spectra in the two cases of the atomic coordinates 



for YVO3 is now under progress. Second, the intensities 
in the (100) reflection are 1/5 ~ 1/10 smaller than those 
in the other reflections in the experiment of YVO3 (and 
YTiOs). Note that the samples are different for differ- 
ent reflections in the experiment. This may make the 
comparisoQ of the data between different reflections less 
accurate.E3 In any event, this behavior is not reproduced 
by the present calculation. It might improve the situa- 
tion to go beyond the muffin-tin approximation in which 
the potential is crudely treated as a constant in the in- 
terstitial region. Third, in addition to the RXS spectra 
mentioned above, another small peak has been observed 
about 10 eV below the main K edge peak. Since its az- 
imuthaJ angle dependence is nearly the same as in other 
peaksjj this pre-edge peak must come from the dipolar 
process. The presence |Qf such pre-edge peaks has been 
predicted for LaMnOsO This study suggests that the 
pre-edge peak arises from the p symmetric states (cen- 
tered at a Ti site) hybridizing to neighboring 3d states 
in the intermediate states. Since the orbital polarization 
is obtained very small in the present calculation, we are 
unable to handle the pre K edge structure. The LDA-I-J7 
method may be suitable to treat this problem. 
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